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We study the dynamical stripe instabilities on the geometries with hyperscaling violation in the IR, 
which asymptotically approach AdS4 in the UV. The instabilities break the translational invariance 
spontaneously and are induced by the axion term ~ aF A F in the bulk action. We first study 
the perturbation equations in the probe limit, and find that there is a strong correlation between 
the stripe instabilities caused by the axion term and parameters of the theories which determine 
the IR hyperscaling violation. Contrary to the IR AdS2 case, the effect of the axion term for the 
stripe instabilities can be enhanced/suppressed at low temperature depending on the parameters. 
For a certain one-parameter family of the hyperscaling violation, we find the onset of the stripe 
instability analytically in the axion coupling tuned model. For more generic parameter range of 
hyperscaling violation, we study the instability onset by searching for the zero mode numerically on 
the full geometries. We also argue that quite analogous results hold, after taking into account the 
graviton fluctuation, i.e., beyond the probe limit. 
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I. INTRODUCTION 

Holography is an extremely useful tool to understand 
the quantum field theory in the strongly coupled limit [H- 
13] . This motivates us to apply the holographic technique 
to interesting subjects such as QCD and condensed mat- 
ter physics, and see what we can learn about real world 
physics from holography. 

In studying the quantum field theories for QCD or con- 
densed matter physics, one of the most important things 
is understanding the infrared (or vacuum) structure. For 
example, in many of the recent studies of the non-Fermi 
liquid in holographic setting, the most interesting nature 
is originated from its IR behavior, i.e., the bulk black 
brane near horizon geometries [IHS] ■ This is natural from 
the trivial fact that Fermi surface is the low energy na- 
ture of the systems. In the holographic superconduc- 
tor case [7HT^. with or without various lattice effects 
|13H15j , the most interesting features are their symmetry 
breaking pattern in the IR, i.e., in the near horizons of 
hairy black branes. Another interesting aspects of black 
brane geometries for the condensed matter application 
is the existence of flux at the horizon, since this gives 
the violation of the Luttinger theorem for the field the- 
ory duals [T6H2T] and this resembles the fractionalized 
Fermi-liquids [1^1 HH 113] . Therefore it gives the possibly 
interesting dual of them. Quantum criticality is another 
fine example where IR dominates the physics. In these 
ways, we are especially interested in the IR behavior of 
the black brane near horizon geometries in holographic 
condensed matter, and it is quite interesting to under- 
stand their various dynamics from the geometry side. See 
also fMI, for recent development of the classification 
of the black brane geometries which admits homogeneous 
but anisotropic geometries. 

One of the most interesting geometries in this view 
point which people actively studied recently is Lifshitz 
geometry [26H28] and so-called geometry with hyperscal- 
ing violation [55H3H]- These geometries are interesting 
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since they are, if realized in the IR, dual to the field 
theories which break the Lorentz-invariance but respect 
spatial rotation and translational invariance in the IR. IR 
Lifshitz and hyperscaling violating geometries can be re- 
alized as the near horizon geometries of some black brane 
solutions, and they can be obtained, for example, in the- 
ories where dilaton has run-away behavior governed by 
the exponential potentials [28-30,- 

Another interesting nature of these geometries is that 
these geometries do not admit large entropy unlike 
Reissner-Nordstrom black branes |55| . On the boundary 
field theory side, this means that it has vanishing entropy 
density at the zero temperature limit, which is more nat- 
ural from thermodynamical view point. In addition, rich 
behavior of the hyperscaling violating geometries allows 
more exotic behavior for the fermion correlators. For ex- 
amples, it gives various w-dependence for the non-Fermi 
liquid decay-rates, as shown in f3(F. 

In the meantime, recently a lot of recent progress 
was made in understanding the translational symmetry 
breaking in the holographic setting [40H47] . The symme- 
try breaking can be induced by the axion term aF /\F in 
4d gravity (where a is neutral pseudo-scalar field, and we 
call it axion in this paper), and can occur both sponta- 
neously and by the source term. In many of the sit- 
uations, the analysis is mainly done for the Reissner- 
Nordstrom black brane which admits AdS2 geometries 
in the IR. Given this, it is very natural to study if above 
IR Lifshitz or geometries with hyperscaling violation can 
survive after we take into account such axion effect and 
see if it induces the translational symmetry breaking. 
Furthermore, since this instability changes the IR na- 
ture of the geometries, studying this instability on the 
hyperscaling violating geometries is by itself interesting 
questions as general relativity problem. In this paper, 
we study the instability on the IR hyperscaling violating 
geometries, focusing on the onset of the stripe instability. 

The organization of this paper is as follows; We first 
review the geometries with hyperscaling violation in the 
4d Einstein-Maxwell-dilaton-axion system in §11. The ge- 
ometries with hyperscaling violation emerge in the IR, 
i.e., in the near horizon limit of the full solutions which 
approach AdS4 in the UV asymptotically. In §111, we 
study the perturbation on these geometries in the probe 
limit. We first analyze the perturbation focusing on the 
IR geometries at small nonzero temperature in §111 B[ 
and at zero temperature limit ^ III C 1 These studies al- 



low us to find the dependence of the stripe instabilities 
on the parameters of the theories which determine the IR 
hyperscaling violation. Given this, we identify the insta- 
bility onset on a certain one-parameter family of the hy- 
perscaling violating geometries, in §111 C 2| In §111 D[ we 
search for the zero mode on the full geometries for more 
generic parameters numerically. In §IV, we study the 
perturbation equations without taking the probe limit, 
and see that essentially the same results hold compared 
to the probe limit case. We end with a summary and 
discussion in §V. 



II. EINSTEIN-MAXWELL-DILATON-AXION 
SYSTEM 

A. The set-up 

The action we consider is Einstein- Maxwell theory cou- 
pled to a dilaton-axion, given by 



S ^ I d^x^^g{R - 2(V0)2 - 2e2«'^(Va)2 - V{(t)) 

-f{q^)F^,F'^'^ - e{a)F^,F^A . (1) 



Here, Ff"" = ^e'^'P'^Fp^, and e^"'"' has a factor of 1/^^ 
in its definition such that axion term 6{a)Fp_^FP-'^ is in- 
dependent of the metric. We take the convention that 

^trxy ^ CI- 
As an explicit example, in this paper we consider 

= e^""* , = 2Focosh2(5<?!), 9{a):^cia, (2) 

for the explicit stability/symmetry breaking by taking 
various real parameters a, (5, ci, and Vq, but we take 
Vq < 0. Note that ¥{(/)) Voe^^'*' in the J ^- oo limit. 
The Einstein equation in trace reversed form becomes 

= (^F^xF^^ - Ig^^.F^^ + l9,..V{^,a) , (3) 

which is irrelevant to the axionic 9{a)F A F term. The 
equations of motion for dilaton and axion are 



--d^i^gg'^^d^a) = {dae{a))FF , 



(5) 



(6) 



(7) 



and for the gauge field, we have 

with the Bianchi identity 

d^,F,x + d.Fx^ + dxF^, = . 

B. Background solutions 

1. IR Hyperscaling violating geometries 



First we review the background hyperscaling violat- 
ing geometries [29l [30] , on which wc later add the small 
fluctuation to study its stability. We consider the back- 
ground geometries where all quantities are functions of 
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only the radial direction r. By using the r coordinate 
re- definition, we make the metric in the form 

ds^ = ~&{rfdt^ + ^ + b{rf{dx^ + dy^) (8) 

and equations of motion and the Bianchi identity are 
solved by 



F ^ Ftrir)dt A dr , Ftr{r) 



Qe 



(9) 



In the IR limit, we choose 6 cj) ^ oo, such that V{(p) — > 
^g250_ Then with the approximate potential V{(f>) — 
yj^g250 g^j-^ obtain the solutions as given by 

EQl, 



-airf = ( 1 

b{r) 



^2 _ ^ 



r / 

r) = k log r . 



where 



k = - 



4 + (a + (5)2 
2{a + 6) 



7 = 1 



4+ (a + (5) 



2 ' 



Qi = -Vc 



2dia + S) 
4+(a + 5)2 ' 
2-5(a + (5) 



2(2 + a(a + 5))^ 



(4 + (a + 5)2)^ 



2 (2 + a{a + 5)) (4 + (3a ~6){a + S)) 



(10) 
(11) 

(12) 
(13) 

,(14) 



and axion field a is set to be zero, a — [56 . There are 
parameter space which satisfies > 0, > 0, 7 > 0, 
7 — /3>0, /3>0 and Sk < 0. The condition 7 > arises 
from the requirement that gtt vanishes at the horizon 
in zero temperature limit, and 7 — /3 > arises from 
the null energy condition at the rji — case. /3 > is 
automatically satisfied and S k < is for (5 ^> 00 at 
r — )■ 0, which gives S{a + S) > 0. 

The temperature T of the system is given in terms of 
rh as 



T 



(27 + 2/3- 1)C/ 27-1 
47r '^'^ ■ 



(15) 



So we require 7 > 1/2 so that at T — > 0. Then, 

with /3 > 0, this geometry is thermodynamically stable 
since the horizon area vanishes at the zero temperature 
limit and dual theories satisfy the 3rd law of the thermo- 
dynamics. 

In zero temperature limit, the metric reduces to 



a(r) = Car^ , b{r) = , 



(16) 



therefore the horizon locates at r = 0. This zero temper- 
ature metric is sometimes written as 

dsl = (-f-2(--l)df2 + ;2^^2 ^ ^^2 ^ ^^2^ ^ 

7-/3 n A ( '^^l 



z = l- 



1-7-/3 
1 



= 4 



1-7-^ 



(l-7-/?)2C2' 



(17) 
(18) 



by the coordinate transformation f = ^, where d = 

4. 

These are the so-called "geometries with hyperscaling 
violation" . The parameter z is called dynamical critical 
exponent and 9 is called hyperscaling violation parame- 
ter. For generic values of a and 6, we have 7^1 and 
6* 7^ 0. If 7 = 1, and equivalently, 9 = 0, this IR met- 
ric has additional scale invariance symmetry, therefore 
9 characterizes the "deviation" from the scale invariant 
limit. Similarly, if z = 1, then time and spatial coordi- 
nate transform equivalently, therefore z characterizes the 
"deviation" from the relativistic limit. 

Before we continue, let us pose to check the number 
of the parameters of our system. Our model, given by 
eq. ([T]) and has 5 real parameters, a, S, Vq, ^, ci. 
However Vq is a parameter to set the scale (of asymptotic 
AdS4), so we can set Vq = —1, without loss of generality. 
Furthermore, since our background has zero VEV for the 
axion field, a = 0, parameters ^ and ci are irrelevant for 
the background. In this way, the background is charac- 
terized by two real parameters a and S only. Especially in 
the IR, these two parameters determine the hyperscaling 
violation 9 and z, or equivalently /3 and 7. 

We will now review the full solutions which interpolate 
these IR /3, 7 geometries with hyperscaling violation to 
UV AdS4 asymptotically. 



2. Full solutions interpolating IR hyperscaling violating 
geometries to UV AdS4 

The solution we reviewed is the IR limit of the full 
solution which asymptotically approaches AdS4 in UV. In 
order to obtain the full solution which connects the IR to 
UV, we perturb the IR solution by 0(e) and numerically 
follow the evolution of it. To obtain asymptotic AdS4 
with appropriate boundary conditions, we tune e, where 
e is a small parameter. In the zero temperature case, we 
set [30] 



a(r) = Car'' (1 + edir-") , 
b{r) = r^{l + ed2r'') , 
0(r) = k log r + ed^r^ , 



(19) 
(20) 
(21) 



and the perturbation equations of motion are all satisfied 
up to the 0(e) by setting 



Vl=V2^V, 0,2 = —di , 
-D2 

, _ A{-1 + v) + [a + 5f{l + v) , 
4(q; + 0) 
2((5(a + (5) + 2) 3 A 
^2^ 



A 



(a + (5)2 + 4 2 2 ((a + 5)2 +4)' 

2 



(22) 
(23) 
(24) 



= ([4 + (3a - (5)(a + ,5)] [(a + 5f + 4] ' 
X (^36 - (a + (5)(a(8(5(a + 5)- 19) + 17(5)^ ^ ^'"^ 



(25) 
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s[r](red) & 
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10 r 




FIG. 1: The interpolating solution \/—gtt{r) = y/g^^ir) = 
a(r) (red) and \/ gxxir) = \/ gyy(r) = b{r) (yellow) for the 
parameter choice a — —1/3, 5 — 0.55, di — 0.01, Vq = —1. 
It shows linir^oo <i{r) — > ^^'^ linir-^oo b{r) oc r as is 

expected from asymptotic AdS4 in the UV. (color online) 



o'[rl(red) & 
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FIG. 2: The interpolating solution 2'(r) (red) and fe'(r) (yel- 
low) for the same parameter choice. limr-»oo a'(T') — > 
(color online) 



1.0 r 



Bi = (a + S) (5a^ + 'Sa^S - 9ad^ + 32a - 76'^ - 16(5) 
+48+ ([4+ {3a- S) {a + 6)] [{a + S)^ + 4^ 

X (36 - (a + S){a{8S{a + 5)- 19) + 17(5))y^^ (26) 
^2 = 2(2 (a^ + 2)6^ + a {3a^ + 4)S-A{a^ + 2) 

-aS^^ X {{a + 5)^ + 4) . (27) 



With this data, we can obtain the interpolating solutions 
which smoothly connect the IR hyperscaling violating ge- 
ometries to the UV AdS4 geometries. To connect our 
analytic IR hyperscaling violating solutions to an asymp- 
totically AdS4 metric, we numerically solve the Einstein 
equations from r = tq to r = 00. Here is some small 



radius where the solutions are given by eq. ( 19 ) and ( 20 ). 
We solve numerically by setting Vq = — 1. At r = tq, the 
condition je^.^rQ | ^ 1 (i = 1, 2) needs to be satisfied [57] . 
We can obtain the appropriate boundary condition in the 
UV by fine-tuning e. 

The results of interpolating functions a{r), b(r), 4>{r), 
for the parameter choice a = —1/3, S = 0.55, edi = 0.01 
are shown in Fig. 1-3. For this parameter choice, we 
obtain 7 « 1 - 0.06, /? « 0.01, k « -0.1, Ca ^ I, Ql ^ 
0.5, and we have positive ly, v 0.9 > 0, and real d2 and 
^3, d2 ~ — 1.4di, da « 0.5 di. The results shown in Fig. 1 
- 3 are essentially the same as given in the appendix F of 



Similarly in the finite temperature case, we can obtain 
the interpolating solutions. We now add the x-dependent 
perturbation on this background to study the stripe sta- 
bility on the background geometry. 



FIG. 3: The interpolating solution 0(r) for the same parame- 
ter choice, lim, — ^ao 0(7") — >■ 0, for the asymptotic AdS4 in the 
UV. (color online) 



III. PERTURBATION ANALYSIS IN THE 
PROBE LIMIT 



A. Perturbation equations in the probe limit 

Given the full solutions which interpolate between the 
IR hyperscaling violating geometries and UV AdS4, we 
add the x-dependent perturbation and study the onset 
of the stripe instability. First, for the sake of simplicity, 
we neglect the graviton fluctuation, namely, we consider 
the stability in the "probe limit" , where matter fluctu- 
ation will not back-react to the graviton fluctuation. In 
this limit, the translational symmetry breaking can occur 
and there can exist unstable mode, but the background 
geometries are fixed. We will see even in this case, we 
have interesting phenomena for a dynamical instability 
of the translationally invariant vacuum. 

On the background metric and flux given by eq. (|8|) 
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and eq. (|9|, we add the following small fluctuation, 



SA, 



5a 



(28) 



and we choose these fluctuations are x-coordinate depen- 
dent. If any of these modes show nonzero vacuum expec- 
tation value with the boundary condition that their non- 
normalizable modes disappear at the boundary, we have 
gravity dual of the translational symmetry breaking. 
From the dilaton fluctuation equation eq. Q, we have 



HdlfmFHcj, -f 2{d^f{cl,))F^^^5F^, . (29) 

Note that the third term on the RHS vanish, since F**" ^ 
only for F*'' but we excite only 5 Ay . Therefore, dilaton 
equation of motion can be trivially satisfied by setting 



= 0, 



(30) 



so we can set the dilaton fluctuation to be zero. In this 
paper, we will set the fluctuation of (j) to be zero, 50 = 0. 
In such case, the operators which are dual to Ay^ a can 
condensate. 

The axion fluctuation equation of motion eq. ([s]) gives 



'-9 



(31) 



Here we have used the fact that background flux is purely 
electrical, FF — 0. 

The gauge field equation of motion eq. ^ yields, 



d^[^-g{fmF^^'' + e{a)5F^^^ 
+{dae{a))F'"'Sa 



0, 



(32) 



here we used S(f> = 0. We need to study the coupled 
equation for SAy and 6a given by eq. (31 ) and eq. (32). 
We follow the analysis of [HI closely. We set, 



SAy = 6 Ay {r) sin qxe ''^ 
5a = (5a(r) cos qxe~*"* . 



(33) 
(34) 



By using the background metric and flux given previ- 
ously, after a bit algebra, the v = y component of eq. ( 32 ) 
gives, 

V=55^^/(0) {-g^'^e - 5"^') 5Ay{r) 
+dr {^ggyyf{<l>)g'-'dr5Ay[r)) 

-q^/^{dae{a))F''y5a{r) = . (35) 



With the metric ansatz eq. ([8|), and by using 5F^y = 0, 
and 



-gF 



xy 



-5e 



xytr TTi 



62/(0) 



(36) 



this equation is re-written as 



/(0)(-r2g2 + a-2^2) 



+dr {f{^)~a'dr5Ay{r)) 
qQ 



-{daeia))- 



5a{r) = . 



(37) 



'62/(0) 

Note that both a term proportional to g^, and a term 
induced by axion term, have explicit common factor b~^. 

Using e''*^y = -e'^^'v = 1/V=5 and the fact that the 
background axion takes the zero VEV, it is straightfor- 
ward to check that v ^ y component of eq. ( 32 1 is auto- 
matically satisfied. 

Similarly, the axion fluctuation equation of motion is 
written as 



+ a,((e2«'^62)a2a,(5a(r)) 

+(5<^%))^'5A.(r) = 0. 



(38) 



We will solve the coupled equations eq. ( |37[ ) and 
eq. ( 38 ) to study the translational symmetry breaking 



and stability of the translationally invariant vacuum. 



However since eq. ( 37 ) and ( 38 1 are coupled, generically it 



is quite difficult to solve these for the stability. However 
since there are analytical expressions for the metrics in 
the IR, which are the hyperscaling violating geometries 
given by eq. (10) and (11), we first restrict our attention 



to the near horizon IR geometries. 



B. Negative momentum square mode on the 
hyperscaling violating geometry: finite temperature 
analysis 

First, let's consider the finite temperature case, where 
the near horizon metric is given by eq. (10 1 and (11). In 



the near horizon region, where r — ^ r^ , we can approxi- 
mate the metric as 



;;2 



(2/3 + 27-l)C^,f-'(^-^/0 



,2fi 
h ' 



/(0) 



„2ak 



(39) 
(40) 



With this approximation, the eq. ( 37 ) and eq. ( 38 ) reduce 
to 



(-r-^^^ + 5-W)M,(r) 
{rl-^a^dr5Ay[r)) + c,^^6a{r) = 

2^k+2P ( -213 2 , ~-2 2\ r / \ 

V [^"^h q +a w j 5a{r) 
+a,(^2Cfc+2/3~2 



(41) 



(42) 



(r^- ■ ' a dr5a{r)) + ci 



^I^My(r) =0. 
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It is more convenient to re-write these as 
qQe 



-ci- 



,2/3+4Qfc 



6a = 0, 



(43) 



qQe 



+Ci- 



(5A„ = . 



(44) 



By multiplying C2 = rjj 



_ B+^k-ak 



to cq. (44), we obtain 



+^i^&^^2=0, (45) 



+Ci^^^i-0. (46) 



where 



My = '01 



C2 (5a = ■(/'2 



(47) 



By choosing ip± = ipi ± ip2, we have two independent 
equations 



±ci- 



.3/3+3afc+^fc 



^± -0. (48) 



These equations have very similar structure to the equa- 
tions of motion for massive scalars in the finite temper- 
ature black brane. However due to the axion term pro- 
portional to ci, it has an "effective momentum square" 
defined as 



q Tci- 



qQe 



fi+?,ak+ik 



(49) 



We take ciqQ^ > without loss of generality. Note that 
this guarantees that at the horizon, one of the "effective 
momentum square" , qeff+ above for the mode 4'+ ('') j 
becomes always negative. The minimum value is given 
at critical wave number qc, 



9c ~ r,B+3ak+^k 

with minimum momentum square 

4Ql 



4r 



2(l3+3ak+(,k) 



T <0, 



(50) 



(51) 



which is always negative. 

This is very similar to the situations where the striped 
phase instability is studied by [lOl HI] ; The axion term 



ci ^ makes the system have lower energy by having the 
nonzero momentum "stripe" type of the spatially modu- 
lated mode. 

However contrary to the IR AdS2 case, there are differ- 
ence on the study of the instability of the IR hyperscaling 
violating geometries with finite but small temperature; 
Due to the additional warping factor in hyperscaling vi- 
olating geometries, there is extra temperature (r/i) de- 
pendence for this effective momentum square, and CiQf. 



^fS+Sak+^k 



appears only in the combinations of ciQe/' 
In the AdS2 case, where Q! = 5 = Ofor/3 = fc = and 
7 = 1 in eq. (12) and (13), ciQe does not involve any 



temperature dependence, since 7.^+'^"''+?'^ becomes con- 
stant. 

If ciQe/r^^^"'''^^'' are very large, for q ^ 0, ip+ mode 
can have lower energy with nonzero momentum q mode, 
and the existence of such negative momentum square 
mode gives the possibility that such mode induces the 
instability of the translationally invariant vacuum. 

The fact that ciQe appears in the combinations of 
ciQe/j^f "'"^"'^"''^'^ implies that the situation is very dif- 
ferent, depending on /3 + 3ak + £^k is positive or negative. 
If 



l3 + 3ak + ^k>0, 



(52) 



as we lower the temperature {i.e., as we lower r/,,), the 
instability are more likely to occur even for the small 
values of ciQe, due to the enhancement by the factor 
l/r^"^'^"*^'*'^'^. In this case, CiQe/r^'^^"'''^^'' diverges pos- 
itively at — > and qeff+ goes to negative infinity, and 
these suggest that the instabilities are expected to occur. 
On the other hand, if /3 + 3ak + ^k is negative, as we 
lower the temperature, we need to increase CiQe in order 
to keep the same value for qlff+, due to the suppression 
by the factor i/^^+3«'=+«fe. 

However, whether the existence of the negative effec- 
tive momentum square mode, qeff+ < 0, which is eval- 
uated at the horizon, is enough or not for instability is 
another question. This is because in general hyperscaling 
violating geometries, we do not know the critical values of 
the instability. This is in sharp contrast to the AdS2 case. 
In the AdS2 case, partially due to the scale invariance of 
the geometry, we have a sharp bound for the instability, 
i.e., Breitenlohner - Freedman (BF) bound. However for 
geometries with hyperscaling violation, we generically do 
not expect sharp bound for generic radius and therefore 
in order to show that there is an unstable mode, local ar- 
gument is not enough in general and we need to study the 
full bulk system to search for the unstable mode. Note 
that since hyperscaling violating geometries are IR limit 
and they approach AdS4 in the UV, this means that we 
need to study the perturbation equations on the whole 
geometries, with the normalizable boundary condition in 
the UV AdS4 region. 

However, before we study the full bulk system for 
the unstable mode numerically, one might wonder if 
something special occurs for the special parameter range 
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where the equahty of eq. (52) holds. To understand this 



marginal parameter range in more detail, we continue 
the study of perturbation equations on the IR hyperscal- 
ing violating geometries by taking the zero temperature 
limit. 



C. Zero temperature analysis and instability 
criteria 

1. Zero temperature perturbation 

We analyze the instabilities at zero temperature limit, 
where the metric is 



(53) 



In this case, by redefining the field as 

5Ay = r^-'-°'^5Ay , 6a = r^-^-''-^^6a, (54) 

the gauge field and axion equations of motion are written 
as 



where 



sx= {fy),M-^ 

Mr) = C^(-l + 7 + afc)(7 + afc) 



Cldydr . 



A{r) B{r) 
B{r) C(r) 



B{r) = -ClgQer•2""'"'^-'^'"+^^ 

Cir) = CU-l + P + l + H)W + l + kO 



(55) 

(56) 
(57) 

(58) 
(59) 

(60) 



where we have set uj 
instability onset. 
Note that if both 



0, since it allows us to see the 



and 



2-3/3-27-fc(3a + < 0, (61) 



2 - 3^ - 27 - fc(3a + < 2 - 2/3 - 27 , (62) 



are satisfied, the off-diagonal terms, which are propor- 
tional to ci and which come from the axion term, dom- 
inate at the extremal horizon, r — 0. In such case, we 
are forced to have the situation where the matrix has 
at least one negative eigenvalue, which goes to negative 
infinity at r — )■ 0. Therefore in this parameter range, it 
is expected that there is an unstable mode for the stripe 
instability. 

On the other hand, if the parameters of hyperscaling 
violation do not satisfy above inequality (61) and (62), 
then the matrix does not always have negative eigen- 
value at r — > 0. However even in this case, if we take 



large ci, then at finite r, has one eigenvalue which 
becomes negative at some intermediate radius r. From 
these, we expect that if inequality ([6T|) and ([62|) hold, the 



instability is more likely to occur compared to the case 
where nequality (61) and (62) do not hold. We will see 



this actually later in the numerical analysis. 

Note that the second condition eq. (62) is equivalent 

to 



l3 + 3ak + £_k>0, 



(63) 



which is the same as eq. (52 1, i.e., the second term of ef 



fective negative momentum square in eq. ( 49 ) , dominates 
at the extremal limit rh — > 0. The first condition eq. (61 ) 



is to guarantee that the off-diagonal term will not vanish 
at the horizon. 



2. Analytical criteria for instability onset for special type of 
hyperscaling violation 

There is a special parameter range in hyperscaling vi- 
olating geometries, which occurs when the equality of 
both eq. (61) and (62) are satisfied. This is when ^, the 
parameter for the axion kinetic term in eq. ([I]), is tuned 
as 



C = -3a - 



3/3 -f 27 - 2 



(64) 



and /3 and 7 in the IR hyperscaling violating geometries 
satisfy the relation 



/3 + 7-1 



(65) 



For this special type of hyperscaling violating geometries, 
we can identify the critical mass just as Breitenlohner 
- Freedman (BF) bound, since the matrix Af^, given in 



eq. (57), becomes independent on the radius r, and there- 
fore constant matrix. From the solution eq. ( 12 ) and (flSl, 



we can see that if we do not have parameter ^, the only 
way to satisfy eq. (64) and (65 ) are AdS2 case, i.e., /3 = 



and 7 = 1 case. However for the theories with tuned ax- 
ion kinetic terms ^ as eq. (64), becomes constant 



matrix for one real parameter family of hyperscaling vi- 
olating geometries with eq. (65), and AdS2 is a special 



parameter point in this parameter range. 

For such IR geometries with hyperscaling violation sat- 



isfying (65), M has two eigenvalues, m 



±1 



where 



\ (^Ao + Co± ^(^o-Co)2+4i32^ 



(66) 



Ao = Ci {-1 + + ak) {j + ak) + q\ (67) 
Bo = -ciqQe, (68) 

Co = cii-1 + p +j + ko i/3 +j + ko + - m 



dX admits the power law ansatz 
SX^r^, 



(70) 



where A satisfies 



C,2A(A + l)=™2±(g), 



(71) 



and critical values for the instability onset are given when 
A becomes complex |59j . By minimizing rn^ {q) with 
respect to q at q = qmin, it becomes 



C2A(A + l)-m2_(g„ 

where qmin satisfies dqm?i{q)\q^q^^, 
complex when 



< 



1 



(72) 

0. A becomes 



(73) 



holds. In the large CiQ^ and q case, 



im?_{q) 



q^ - \ciQe\q, 



(74) 



which allows g™„ « |ci(3e|/2, and m^(q„i„) w 
— |ciQeP/4. Therefore in this case, eq. (731 is satisfied 
and we have stripe instability. 



The special relation eq. ( 65 ) is the point where both 



z and 9 diverge [z and 6 are given by eq. (17l). More 
precisely, it is the special point of hyperscaling violation 
where either z — )■ —oo and 6 — > +oo, or z — )■ +cx) and 
9 — > —CXI holds. 

It is interesting to see how above results are different 
from the AdS2 case. Let's first consider the AdS2 case, 
which is the special limit of above and we take a = 5 = Q 
so that we have = 0, 7 = 1, fc = in eq. (12), and 
eq. (13). In such case, Aq, Bo, Cq in eq. (67) - (69) 



simplifies and we obtain 



"^±AdS2(9) = q^±ciqQe 



(75) 



This iTi^j^^g2 is the same effective momentum square 
leff^ given in eq. (49). Therefore in AdS2 case, the 
onset value of the instability, which we call Cmin AdS2 1 is 
given by the equality of eq. ( 73 ) as 



^min AdS2 



Ca\ 



Qe\a=S=0 



(76) 



Here we have used eq. (13) and (14). The stripe instabil 



ity occurs at Ci > Cmin AdS2 ■ 



Next, we consider non-AdS2 case. Note that eq. (65) is 



satisfied when a = ±S from eq. (12) and (13). However 



since a = —S gives /3 = 0, 7 = 1, which is AdS2, we 
consider a — 6 case. From eq. (12) - (14), in this case 



the hyperscaling violating geometries are parametrized 
by only one real parameter a, with 
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a 



1 



= -Vn 



2(1 + a2) 



k 



Cl = -Vo ., 



(77) 
(78) 



Then, the eigenvalue m?. (q) of given by eq. ( 66 ) be- 
comes 

2 (a'^ + 1) ^ 

-V2^J (a2 + if Vo ((a* - 1) cfq^ + Sa^Vo)) . 

(80) 



Since this has minima at 



2\2 2 _ 64"! 



2^2(1 -a4) 
the instability condition eq. ([73| becomes 



(81) 



8(a2 + l)' 

(a2 + l)Vo 



- x((a^-l)^(a^ + l). 



8c? (a2 - 1) (a2 + 1) 
-64 (a^ + a^) - 32 (a^ - l) a^cl) < 



Vo 



(82) 



Therefore, the onset ci value of the instability, which we 
call c^j„, is given by the equality of this and it is 



V2 



1 + a' 
1-a^ 



2\ 1/2 



(83) 



The instability occurs when ci > c^j„, which is larger 
than the critical value CminAdS2 AdS2 case, 
Therefore it is more stable than the AdS2 case, regarding 
the stripe instability. 

This expression loses its meaning when > 1, cor- 
responding to 7 < 1/2. But 7 cannot be less than 1/2, 
since then the entropy density of these black brane are 
proportional to the negative power of the temperature, 
which, thermodynamically, does not make sense |60j . An- 
other reason for 7 > 1/2 is that if 7 < 1/2, then the null 
rays from the zero temperature horizon r = can reach 
nonzero r at finite time, which contradicts with the "hori- 
zon" property at r = [5T] . 



3. C = case 



In the next section, we conduct numerical analysis of 



generic parameter points where such an equality eq. ( 65 ) 



does not hold. For that purpose, let's investigate the 
^ = case a bit more, since this is the case where axion 
has canonical kinetic term. 



by 



The parameter range which satisfies eq. ( 62 ) is given 



{a + 6) (5a-5) <0, 



(84) 



with Vo < 0, and eq. (64) gives 



and for eq. (61 ), it is given by 



^ = -2a 



(79) 



(a -I- S) (3a + S) <0 



(85) 
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This gives 



3a < 5 < —a (for a > 0) , 
-3a > S > -a (for a < 0) . 



(86) 
(87) 



This range is written in terms of /3 and 7, or 9 and z. 
From eq. (12 1, we have 



^ 1-2/3-7 . , 7-1 
a — ± — , , = ±- 



V/3(l-/3) ' " V/3(l-/5) 
Remember that we need 

^ < 7, < /3 < 1. 



(89) 



Therefore, the parameter range, eq. (86) or (87 1, is equiv- 
alent to 



0</3<^(l-7), ^<7<1. 



In terms of 9 and z, using 



9-4 _ 0-Az 

^ ~ 2{9-2z) ' ^' ^ 2(6*- 2z) 



(90) 



(91) 



from eq. (17), the parameter range eq. (901 is written as 



4<6i<6, z<0. 



(92) 



So far we consider the parameter range satisfying 
eq. (61) and (62). In addition to these, we have more 



restriction on the parameters. First, we need 5k < Q 
which restricts 5{a + (5) > 0. Second, we want to study 
the stripe instabihty on the background which is stable 
at (J = 0. Since the IR hyperscaling violating geometries 
asymptotically approach AdS4 in the UV, the mass of 
the dilaton must be above the AdS4 BF bound for the 
stability at zero momentum q = Q. From our Lagrangian 
eq. ([1]) , we can read off the mass of the dilaton at = 
and this gives additional constraint 



<5< 



(93) 



Combined these additional conditions with eq. (86) 
and ( 87 ) , finally they are summarized as 



> -]^5 > a> 



-5 > 



or, 



(94) 



(95) 



It is clear that there are parameter ranges satis- 
fying above in the "physical parameter ranges" [62] 
parametrized by a and 5, as is shown in Figure 1 of |30j . 
In this parameter range, the effective momentum square 



in eq. (51) becomes large negative value at low tempera- 
ture, and also the off-diagonal components of matrix 



in eq. (57) dominates. Therefore it is expected that the 
instability is more likely to occur at small value of ci in 
this case. 

We next conduct numerical analysis for the parameter 
choice where eq. (94) is satisfied and also for another 
parameter choice where it is not satisfied. We call the 



case where eq. (94) is satisfied as case I, and the case 



where eq. ( 94 ) (or eq. ( 63 ) in more generic case where 



^ ^ 0) is not satisfied as case II. 

The limit a — > and J — ?> corresponds to the AdS2 
and its stripe instability is well studied in [HI • We 
will now study numerically the whole bulk system in 
these parameter range for the stripe instability next. 



D. Numerical analysis for the bulk zero mode on 
the full geometries 

Now we conduct numerical investigation of the coupled 
equations eq. (37) and eq. (38) to study the dynamical 
translational symmetry breaking. So far we have been 
concentrating on the perturbation equation analysis re- 
stricting our attention on the IR hyperscaling violating 
geometries. However, now we will conduct the numerical 
analysis on the full geometries, which approach geome- 
tries with hyperscaling violation in the IR and AdS4 in 
the UV. 

In order to find the onset of the instability, we will look 
for the zero mode, namely the solution of eq. (37) and 



eq. (38) with w = 0. If there is such a zero frequency 



mode for some value of Ci, then it is the critical mode. 
By increasing ci above the critical value, the instability 
occurs. This is because if there is an unstable mode, 
Im((u;) > which grows as time evolves, then there should 
also be zero mode solution a; = at the instability onset 
point. 

Let us emphasize again why searching the zero mode 
on the full background geometries is important: In the 
case of IR AdS2, we have a sharp local criteria for the 
onset of the stripe instabilities, which is given by the 
condition ci > Cmm AdS^ with eq. ( 76 ) . This is in contrast 



to our generic hyperscaling violating geometries in the 
IR; in two-parameter (/3 and 7, or 9 and z) hyperscaling 
violating IR geometries, we do not have sharp criteria for 
the stripe instability from the IR geometries. 

However even in the IR AdS2 case, it is important to 
find the zero mode on the full geometries which approach 
AdS4 in the UV. This is because generically the stability 
is not determined by the IR region only, but rather it is 
determined by the full geometries. It is especially so if 
two modes are coupled. 

To see this, let us consider the IR AdS2 geometries 
which approach AdS4 in the UV. For the IR AdS2 region, 
the matrix in eq. ( 55 ) becomes constant matrix and 



we can obtain two eigenvector modes (let us call these as 
mode A and B) made by some specific linear combination 
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of mode SAy and Sa. Suppose that the mode A gives 
the lower eigenvalue of the matrix than the mode 
B, and furthermore that the eigenvalue of the mode A 
is lower than the AdS2 BF bound and the eigenvalue 
of the mode B is higher than the AdS2 BF bound. It 
is true that if we excite only the mode A, then wc can 
lower the energy by the mode A perturbation and this 
indicates instability. However due to the UV boundary 
condition and the mixing of the two modes, generically 
we cannot excite only the mode A, but rather we need 
to excite the mode B too in general [S3]. At what ratio, 
we excite the mode A and B, is determined by the full 
(whole) geometries. Therefore, depending on the ratio of 
the mode A and B excitations, we can see if the system 
is unstable or not by the perturbation. This is because 
the mode A lowers the energy but the mode B increases 
the energy in the IR AdS2. In this way, it is generically 
determined not by the local IR geometry only, but rather 
by the full geometries. Therefore it is important to find 
the zero mode on the whole geometries including IR and 
UV. 

Furthermore, in the case of IR hyperscaling violating 
geometries, which is parametrized by two parameters, we 
do not generically have sharp BF-like bound. Therefore 
it is more important to search for the zero mode on the 
full geometries which include both IR and UV region. 

We now seek for the zero mode on the full geometries. 
For numerics, we set Vq = —1. Introducing new variables 
-0-1- = SAy±C2ip2 as eq. (47), we can obtain the boundary 
condition at the horizon r = hy imposing regularity 
as 



drtp 



±\r= 



1-2/J-27 2 „/. I 

(2/3 + 27-l)C2 



(96) 



where qlff± is the effective momentum square defined in 
eq. (|49|. 



In order to find the onset of the spontaneous transla- 
tional symmetry breaking in the holographic dual setting, 
we impose Dirichlet boundary condition for the variables 
SAy and Sa, at the AdS4 boundary. 



SAy = Sa = 



for 



(97) 



This corresponds to the requirement that there is no 
non-normalizable mode for SAy and Sa in the UV AdS4 
boundary. Non-normalizable mode for SAy and Sa ap- 
proach constant in the UV AdS4. 

It is useful to consider the parameter counting. For 
given fixed parameters a, (5, ^, r^, and ci, there are two 



and q. Note that V 



+ \r=rh 



free parameters V'-lr 
not a free parameter, since we can fix ip-f-lr^rh to unity 
without loss of generality in linearized perturbations. We 
tune ?A_|r=rh and q such that the two boundary condi- 
tions eq. (97) are satisfied. After this tuning, we have 



no parameter left, and as a results, we have nonzero nor- 
malizable mode for both and i/'- at r — > oo. This 
implies that there are spatially modulated VEV for the 



scalar and vector current < jy >, which are dual to ax- 
ion and and gauge boson Ay, and this implies that dual 
theories at IR show the "current density wave" phase. 
Note also that for a given temperature T(r^), we expect 
normalizable zero modes to appear at specific values of 
momentum q. 

As we have discussed, we expect that the effect of 
the axion term and the negative "effective momentum 
square" is enhanced or suppressed depending on eq. ( 94 ) 
is satisfied or not. In order to see the difference, we 
investigate two typical cases; Case I corresponding to 
a = -0.33, S = 0.55, C = 0, which gives /? = 0.01, 
7 = 0.94, k = -0.10, and it satisfies /3 + 3ak + £_k > 0. 
Case II corresponding to a = 0.2, S — 0.55, ^ = 0.4, 
which gives /? = 0.12, 7 = 0.82, k = —0.33, and it satis- 
fes (3 + 3ak + ^k< 0. 

It is easily checked that for the case I, it satisfies 
eq. (|94|, and therefore, ([6T|) - (|63|). So case I corresponds 



to the case where axion term is expected to be enhanced 
at low temperature due to ciQe/r^ 



l3+3ak+^k 



—T' 00, 



and we expect that axion term induces instability even 
at very small values of ci at low temperature. 



On the other hand, for the case II, it violates eq. (94) 



(more precisely it violates cq. (63) since ^ 7^ 0), but does 
not violate the UV AdS4 BF bound (|93l). Therefore case 



II corresponds to the case where axion term is expected 
to be rather suppressed at low temperature r/j — >• due 
to ciQe/r^'^^"'''^^'' — >■ 0, and we expect that we need 
large values of ci in order to induce instability at low 
temperature. 

For both cases, we first numerically construct the back- 
ground solutions interpolating the IR analytical hyper- 
scaling violating geometries and UV AdS4 as we re- 
viewed in jjll B 2 Then we solve the coupled eq. (37) 



and eq. i38l and find numerically the solutions 4'± satis- 



fying the boundary conditions eq. ( 97 1 for a suitable ci , 
q and ip- for each temperature T {rh), where eq. (15) 
gives the relationship between and T. 

For each temperature we find the zero mode solution by 
tuning q for ci , as far as ci > Cmin ■ There is a minimum 
value Cmin for ci, namely, we could not find any zero 
mode solution for any q and ip- when ci < Cmin- 

We plot Cmin for each temperature T in Fig. 4 (case 
I) and Fig. 5 (case II). As expected, for the case I, Cmin 
decreases significantly as T decreases and seems to ap- 
proach zero. On the other hand, for the case II, Cmtn 
does not change drastically as T decreases. This sug- 
gests that spontaneous translational symmetry breaking 
easily occurs for much lower temperature when eq. ( 94 ) is 



satisfied, although we cannot further evaluate the min- 
imum value Cmin for much lower temperature. This is 
because highly numerical accuracy is required. Similarly 
we face another difficulties at logj^Q T > —3.6. This is 
because in such temperature, is not so small and as 
a result, (f>\r=rh is not large enough. Then, our pertur- 
bative method in §11 B 2| to construct full interpolating 
solutions breaks down. Therefore we conduct numeri- 
cal analysis in rather restricted low temperature range. 
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Cmin 

13 - 



• • • 



Logr 



• - 



FIG. 4: The plot of Cmin for various temperatures T for case 
I (a = —0.33, S = 0.55, f — 0). Each plot corresponds to 
(c™„, T, q) = (1.3, 2.05*10"'', 0.83), (1.1, 2.69*10-^ 0.85), 
(0.9, 2.07 * 10~^ 0.86), (0.8, 6.08 * 10"^, 1.0). The parame- 
ters are approximately on the line CminiT) — 2.0+0.20 logjQ T 
in the temperature range we study. 



Cmin 

6 - 

5 - 



- Logr 

-3.6 ° 



FIG. 5: The plot of Cmin for various temperatures T for case 
II (q — 0.2, 5 — 0.55, £, = 0.4). Each plot corresponds to 
(c™„, T, q) = (3.2, 2.36*10-^ 0.45), (3.2, 5.43*10-^ 0.45), 
(3.2, 1.25*10"^, 0.45). c™„(r) is almost constant, c„i„(r) x 
3.2 in the temperature range we study. 



iio ^ —3.6 for case I, and —4.9 < log^o 



T < 



—6.3 < log^ 
—3.6 for case II. 

Note that the value range of Cmin for case I in Fig. 4 
is lower than the critical value for the AdS^ case, 
CminAdS2 = V^, which wc obtained in eq. (76). This re- 



sult is consistent with the analysis in ^ IIIB and ^ III C 1 
that these are instabilities triggered by the enhanced 
axion term effect due to the radius dependent factor, 
ciQe/i'h^^'^^^^'^ — > oo at r/i — > 0. One might expect that 
in case I, at the zero temperature limit where — > 0, 
the critical value Cmin approaches zero. Fig. 4 is consis- 
tent with this expectation. However due to the numeri- 
cal difficulties, we could not confirm this at the very low 
temperature T < 6 * 10^''. 

On the other hand, in case II, we find that the critical 
value Cmm is almost constant as we lower the tempera- 



FIG. 6: T — q curve for the normalizable zero mode in full 
geometries for the case I parameter choice (a = —0.33, 5 = 
0.55, 5 = 0), with Ci = 1.3. Due to the "enhancement factor" 



1/^ 



'h 



8, even though ci is smaller than the critical 



value for IR AdS2 case, CminAdS2 ~ V^, we have two g's 
for the zero model, and in between, there should be unstable 
modes. This figure is in good comparison with IR AdS2 case 
analysis done in 41 , Figure 2, where, ci > CminAdS2- 



ture Th- Note that the value range of Cmin for case II 
in Fig. 5 is higher than the critical value for the AdS2 
case, suggesting that these are instabilities triggered by 
the suppressed axion term effect due to the radius de- 
pendent factor in case II. However we do not have clear 
physical interpretation of the result in the case II. From 
the IR analysis, it might suggest that the critical value 
Cmin increases as we lower the temperature, but the result 
of Fig. 5 is not so. One possible reason for this behavior 
is that the negative momentum square is more dominat- 
ing away from the horizon, r ^ r^- As we have discussed 
in ^ III C 1 one can see that even in the case II with zero 
temperature limit, one of the eigenvalue of can be- 



come negative at some finite radius r. In other words, one 
of the eigenvalue of can become negative at some fi- 
nite radius r, but as r — )■ 0, that value approaches zero. 
Because of this, even though one of the eigenvalue of 
becomes zero at the horizon in the case II, it can become 
some negative value at some finite radius r and therefore, 
there could exist a zero mode in the whole bulk for the 
stripe instability. In such case, if the bulk region, where 
eigenvalue becomes negative, is away from the hori- 
zon r^, then it is possible that lowering the temperature 
does not influence these bulk region much. As a result, 
in such a case, changing the r^ does not change the Cmm- 
However in order to obtain clear physical understanding 
of these results, we need more detail analysis. 

It is also useful to draw the figure for the zero mode in 
{q,T) plot with fixed ci value. [SJ. We plot the critical 
temperatures T versus q for the normalizable zero mode 
in Fig. 6 (case I parameter choice with ci = 1.3) and 
Fig. 7 (case II parameter choice with ci = 3.3). 

In the case I, for given temperature, we generically 
have two g's allowing the normalizable zero modes. Un- 
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FIG. 7: T ~ q curve for the normalizable zero mode for the 
case II parameter choice (a — 0.2, 5 — 0.55, ^ = 4), with 
Ci — 3.3. 



stable modes should exist for momentum q between 
gmin(r) and qmax{T), and at T « 2.05 * lO"'^, Qmin and 
Qmax coincides at g « 0.83. As we lower the temperature, 
<lmin{T) decreases and qmax(T) increases. In the case I, 
because CiQe/r^'^^"'''^^^ goes to infinity at — >■ 0, we 
expect to have very large effect of the axion term, and 
therefore we expect that Qmax becomes very large, at zero 
temperature limit. It would be really nice to confirm this 
numerically, however we could not conduct numerics at 
this very low temperature, T < 4.2 * 10~^ due to the 
difRculties of numerical analysis. 

Again in the case II, we do not have clear physical 
interpretation of the results. 



IV. PERTURBATION ANALYSIS BEYOND 
THE PROBE LIMIT 

We have analyzed so far the system without the gravi- 
ton fluctuation, namely in the probe limit. However, it 
is pretty straightforward to conduct the similar analysis 
with the graviton fluctuation, and we will see that the 
analysis with graviton fluctuation shows essentially the 
same results, compared with the analysis without gravi- 
ton. We will see here that there is a negative momen- 
tum squared mode in the hyperscaling violating geome- 
try even after we take into account the coupling to the 
graviton. 

We consider again the action ([T]). The fluctuations we 
consider are the following components, 6gf 
and we take the following mode dependence, 



'jty, SAy, 



Sgty 

6A. 



y 
Sa 



Sgty{r) sin qx , 
SAy{r) sin qx , 
5a{r) cos qx , 
0. 



(98) 
(99) 
(100) 
(101) 



which has additional graviton mode. Here we have set 
a; = 0, i.e., no time-dependence from the beginning com- 



pared with eq. ( 33 1 and ( 34 1 , in order to discuss the onset 



of the instability 

Quite analogously to the previous analysis in §111 A 



the probe limit, given the background geometry eq. 
we have the equations of motion for the gauge field SAy, 

^f{4>)h-\HAy{r) 

+dJf{cjy)~a^dr5Ay{T) + f{^){drAt)6gty{T] 



-idaOia))- 



-5a{r) = 0, 



(102) 



and for the axion, 

+qida9{a))FtrSAy{r) ^0. (103) 

This equation is unmodified by Sgty ^ 0. 

In addition, from the (t, r) component of the trace re- 
versed Einstein equations, we have fluctuation equation 
for the graviton, 

a^b^ {{d^Sgtyir)) - 4f{cb)Ftr{drSAy(r))) 

+ (4ab{dr&){drb) + b^ {2fmFtrf + V{cl))) 



-q^ygty{r)^0, 



(104) 



and all the other components of Einstein equations and 
equations of motion are automatically satisfied. 

Let's investigate the near horizon limit in similar way 
to the analysis of §111 B| and §111 C 1| From graviton fluc- 
tuation eq. (104), by imposing the regularity condition 



of the solution at the horizon, we can see that we need 



the boundary condition Sg- 



ty\ 



at the hori- 



zon. Then, it is more convenient to set the new variable 
Shty — dr-Sgty And quite analogously to the case where 
we neglect graviton fluctuation in §111 C H the three equa- 
tions eq. (102 1 - (104) are approximated and written in 



the near horizon as 
where 



M^\r.=rJY\r^r 



SY 



SAy 

6a 
Or (Sgty) 



(105) 



(106) 



and 



Af2 = 



<zV62 
(M2i(r))2 



(Mi2(r))2 



Ftr 



(M33(r-))2 



= -{daB{a))qFtr/f{4>), 
= ~ida9{a))qFtr/e'^H\ 

= -A{daOia))q{Ftr)\ 



62 



2,/(0)(Ft,) 



4:ab{dra){drb) 
62 



(107) 

(108) 
(109) 
(110) 

(111) 
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and M^|r=rft means evaluated at the horizon r = rh 
The determinant of AP gives 
DetM^ 

-{da6{a)fq^{Ftrfr\<l^)e-^^*b-\M:i^f 
^ / ,_{dae{a)fqH\Ftrf^ 



xh-'{{A'h:,f~Af{4>){F,rf) . (112) 

From eq. (12 1 and (13), we can obtain the relation 

- 2a/c - 4/3 = 2(5fc = 27 - 2 (113) 

in the near horizon region of geometries r — > r/j with 
hyperscaling violation, where V{<f)) — >■ Vqc'^^'^, then using 
eq. ^ - ([n]), we have 



DetA/2 

= rl'\q^-q^ciQir- 



4 „2 2^2^-2(/3+fc(3a+^)) 



, 2 /• 27+2/3-2 



(114) 



where 



C = 2(Qe)' + 2(2/3 + 27 - 1)/3C2 + Vo . (115) 
Therefore, if 



/3 + fc(3a + > 0, 



(116) 



is satisfied, then the axion term proportional to ciQp 
dominates. This is exactly the same condition, eq. (52) 
and eq 



( 63 1 , which we have obtained in \ III B and 



III CI in the probe limit. In above DetAf^, gauge bo- 
son fluctuation 5Ay and axion fluctuation 5a gives the 
factor proportional to (9* - g2c2Q2^-2(^i+M3a+C))^^ ^^^^ 
that this is q^ffj^ x ^e//- defined by eq. (49). And it is 
essentially DetM^ given in eq. (57) - (60 1, with eq. (61) 
and (62) at r — > limit, up to overall factor r^^^^^^^^. 



Therefore even with the graviton fluctuation, 5 Ay and 5a 
give the same mode for the instability. 

The addition of graviton fluctuation, simply adds one 
more eigenvalue to above Det Af^, and that eigenvalue is 
proportional to [q^ — (rf^^'^^^'^). Therefore, even with 
the graviton fluctuation, the existence criteria of the neg- 
ative eigenvalue mode of the matrix DetM^ due to the 
axion term is precisely the same as the case without gravi- 
ton fluctuation; We expect that if (116 1 is satisfied, the 



effect of axion term are enhanced as we lower the tem- 
perature and intrigues more stripe instability, and the 
behavior of minimum ci for stripe instability is expected 
to show very similar behavior to the Fig. 4. On the other 
hand, if ( 116 ) is not satisfied, we expect that the effect of 



at the horizon and this intrigues less instability, and the 
behavior of minimum ci for stripe instability is expected 
to show very similar behavior to the Fig. 5. It would 
be best if we can confirm this by solving the eq. ( 102 ) 



axion term are suppressed as we lower the temperature 



( |104 1 numerically with the normalizable boundary con- 
dition at the UV AdS4 boundary as we have done in the 
probe limit in §111 D| However in this case, the parame- 
ter range we seek for the normalizable boundary condi- 
tions becomes 3-dimensional, instead of 2-dimensional for 
the probe limit case, and this turns out quite hard task. 
Therefore we leave this as future work on this paper. 

Clearly at the large ci limit, we can have a mode which 
has large negative eigenvalue at some radius, and this in- 
dicates the striped phase instability can occur. One dif- 
ference, compared to the probe limit, is that the analysis 
for the analytic expression for the onset of the instability 
in jj lll C 2 1 does not work. This is because even if both 
eq. ( 64 1 and ( 65 1 hold, the determinant of the matrix 
is proportional to r~^^ . So we need /3 = and 7 = 1 for 
the matrix to become constant matrix, which is AdS2 
case. This implies that we need to introduce one more 
parameter in the Lagrangian to be tuned, so that we can 
have constant matrix in the presence of graviton. 



V. SUMMARY AND DISCUSSION 

In this paper, we studied the stripe instabilities (spa- 
tially modulated instabilities) of the geometries with hy- 
perscaling violation in the IR, which approach AdS4 met- 
ric in the UV asymptotically. The instabilities are in- 
duced by the axion term 5S = J d^xciaF AF in the bulk 
4d action. We first study the perturbation equations in 
the probe limit, and saw that there is a strong correlation 
between the stripe instabilities caused by the axion term 
and parameters of the theories which determine the hy- 
perscaling violation. Contrary to the IR AdS2 case, we 
found that, due to the lack of scale invariance and the 
nontrivial radial dependence of the IR hyperscaling vio- 
lating geometries, the effect of axion term can be either 
enhanced or suppressed depending on the parameters. In 
the parameter range where the effect of the axion term is 
expected to be enhanced, the stripe instability occurs and 
Cmin decreases as we lower the temperature, where Cmin is 
the critical value for the instability and instability occurs 
only at ci > Cmin- On the other hand, in the parameter 
range where the effect of axion term is expected to be 
suppressed, we find that c,„i„ docs not change much as 
we lower the temperature. We have explicitly obtained 
the zero mode solutions for the coupled fluctuations of 
gauge boson 5 Ay and axion 5a numerically in the probe 
limit, with the boundary condition that there are no non- 
normalizable modes. This implies that in the dual field 
theories, the scalar and vector current < jy >, which 
are dual to axion a and gauge boson Ay in the bulk, ac- 
quire the spatially modulated VEV spontaneously, and 
that dual theories at IR show the "current density wave" 
phase. We identify the instability onset on a certain one- 
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parameter family of the hyperscaling violating geometries 
analytically, where the relation eq. ( 65 1 holds. We also 



argue that quite analogous results are expected to hold 
beyond the probe limit. 

There are several open issues which should be under- 
stood in better way. We have done our search of the zero 



mode on rather limited temperature range in i III D This 



is due to the numerical difficulties, and it comes from the 
fact that our background solutions, which are hyperscal- 
ing violating geometries in IR and approach AdS4 in UV, 
are constructed only numerically. If we could construct 
an analytical solution, we can search numerically for the 
zero mode more accurately. So it is interesting and im- 
portant to look for analytical background solutions which 
interpolate between UV AdS4 and IR hyperscaling vio- 
lating geometries. By conducting the numerical analysis 
in better way, we can check if the Cmin goes to zero or 
not in the zero temperature limit in Fig. 4. Similarly, 
we can check how the T — q curve behaves in the zero 
temperature limit in Fig. 6. It is interesting to check 
these. 

We argue in §IV| that the results of instability analy- 
sis are essentially the same by taking into account the 
graviton effects. Of course, it is better if we could con- 
firm this by searching for the zero mode explicitly on the 
full geometries, as we have done in §111 D| in the probe 
limit. 

There are other issues which we would like to under- 
stand in better way. In this paper, we did not argue the 
validity of the action eq. ([T]) for the background solu- 
tions. However, it can happens that the starting action 
is not valid for describing the solutions, depending on the 
behavior of the solutions. For example, our background 
dilaton has run-away behavior, and if the background 
dilaton runs to the strong coupling direction in IR, then 
we have to worry about the possibility that our starting 
action is highly corrected due to the strong coupling ef- 
fects. Such a possibility exists if the starting action is 
derived under the weak coupling approximation. For ex- 
ample, from our action eq. ([T]), the effective coupling of 
the gauge field is given by g'"^^^-^ = f{4>) = e^""^. In 

our background solution this behaves as gu{i) ~^ r""*^ in 
IR. If we assume that our theory eq. ([T]) is derived under 
the weak coupling condition, this with eq. ( 13 1 forces us 



—ak cx a{a+5) > for consistency. But this implies that 
our condition eq. (84) and (85), for the case ^ = 0, can- 



not be satisfied. And we have only the parameter range 
where the axion term is expected to be suppressed, which 



corresponds to case II in the analysis of \ III D These is- 
sues should be understood more from the string theory 
embedding view point. However rather in this paper, we 
study the stability analysis with the assumption that the 
action eq. ([T]) is valid for any solutions, we do not "de- 
rive" our action eq. ([T]) from string theory. It would be 
nice to study these consistency points in more detail. 

It could be that if the background dilaton blows up in 
the IR hyperscaling violating geometries, then we need 
to take into account the higher loop corrections and this 



might make the geometries into the AdS2 metric in fur- 
ther IR, as studied, for examples, in [5TI - f53] . This view- 
point resolves the problematic singularities of the zero 
temperature limit of the hyperscaling violating geome- 
tries at r — ^ [BB]. However, it is not clear if this is 
always the case. For examples, once higher loop correc- 
tions (corresponds to higher string coupling gg correc- 
tions) enter the game, we always need to worry about 
full loop correction effects. Namely, once we face the sit- 
uations where higher order gs effects are as important as 
leading order in gg expansion, this implies that gg expan- 
sion is no more valid. But in general, we do not have such 
a fully non-perturbative effective action, and the validity 
of the loop corrected effective action in order to derive 
the deep IR AdS2 metric is unclear. 

Another interesting question is the end point of the 
stripe instability. In this paper, we studied the onset of 
the stripe instabilities. However to see what is the end 
point of these instabilities, perturbation analysis is not 
enough. We need to study the equations of motion where 
Ay and a are coupled in the probe limit, and the full 
back reacted Einstein equations to go beyond the probe 
limit. For successful examples of the end point of stripe 
instabilities, see gOl EH il] . 

In the probe limit, we identify the onset of the stripe 
instability when the relation eq. (65) holds. It is interest- 



ing to note that this relation also holds at the transition 
point between quasi-particle picture holds/breaks down 
from the study of the fermion Green's function on these 
background [30] ■ It is interesting to investigate to see if 
there are any deep reason for this coincident. 

There are many open questions to be understood 
better. However one thing which is very clear is that 
these geometries with hyperscaling violation and stripe 
instabilities are rich subject and it is worth under- 
standing more in great detail. We hope to return these 
questions in future. 

Note added: When we have almost finished prepar- 
ing for the draft, a paper appeared fS5, where they also 
studied the stripe instability on the hyperscaling violat- 
ing geometries. However our set-up and analysis is differ- 
ent from the one in [SI]. The authors of [SI] studied the 
geometries whose IR (r — >■ 0) is AdS2, and in large r, they 
approach the hyperscaling violating geometries. They 
identified the onset of stripe instability in this IR AdS2 
region by using the AdS2 BF bound, i.e., local AdS2 ar- 
gument. Then they interpret that onset in terms of the 
large r hyperscaling violating parameters. On the other 
hand, in this paper we studied the geometries whose IR 
are hyperscaling violating geometries, which interpolate 
to the AdS4 in the UV. 
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H + D. 



r at r — > 
As shown in 



00, where X± — —3/2 ± 
48j , the theory is generally 
when 
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unstable for the perturbation with respect to 
both modes exist. However for the perturbation anal- 
ysis we conduct in this paper, the dilaton (j> does not 
fluctuate. Therefore we consider that the instability of 
|48| is not related to the stripe instability we investigate 
in this paper. In addition, for the numerical analysis in 
§111 D| we choose D+ as small as possible in the range 
\ed^r^\ < 1 (i = 1,2). 

Qe is actually not a parameter, since it takes definite 
value as eq. ( 13 1. 
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we need to impose boundary condition both at the IR 
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for the second order differential equations, we have im- 
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since Sgty\r=rh — > 0, and 2) using the regularity of 
the flux Ftr and dilaton f{<j}) at the finite tempera- 
ture horizon r — r^, we can neglect term proportional 
to dr5Ay{r) since 

I [f{<t>)Ftr{dra^)dr5Ay(r)] |.=.^ | holds. 



[66] This singularities can be avoided by introducing small 
but non-zero temperature. And our studies of the stripe 
instability in small but non-zero temperature case are not 
affected by this singularities. 



